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Abstract 



We show that the collective vibrational model of atomic nuclei displays 
quantum chaos. To avoid the problems related to the tunneling effects, a 
Green function deterministic numerical method has been used to evaluate 
the energy levels. 
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1. Introduction 

In the last years many studies have been devoted to the properties of 
quantal systems which are chaotic in the classical limit [1] . 

In quantum mechanics, given the impossibility of defining the trajecto- 
ries, classical concepts and methods cannot be applied. Nevertheless, many 
efforts have been made to establish the features of quantum systems which 
reflect the qualitative difference in the behaviour of their classical counter- 
parts. Many schematic models have shown that this difference reveals itself 
in the properties of fluctuations in eigenvalue sequences. The spectral statis- 
tics for the systems with underlying classical chaotic behaviour agree with 
the predictions of the random matrix theory. By contrast, quantum analogs 
of classically integrable systems display the characteristics of Poisson distri- 
bution [2,3]. 

In atomic nuclei, the experimental data of nuclear spectroscopy suggest 
regular states near the yrast line and chaotic states near the neutron emission 
threshold, but the coexistence of regular, chaotic and collective states is not 
yet well understood [4,5]. 

In a previous paper [6] we studied the transition form order to chaos 
in the nuclear roto-vibrational model [7,8], while in this work we analyze 
in greater detail the numerical technique to obtain the energy levels and the 
effects of the order-chaos transition in the spectral statistics for the collective 
vibrational motion. 
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2. The model 



The two-dimensional model which describes the collective vibrational mo- 
tion of atomic nuclei was introduced by Bohr and Mottelson [7] and developed 
by Eisenberg and Greiner [8]. The Hamiltonian is given by 

H=]^B{xl + 2xl) + V{x,,X2). (1) 

where 

y (xi, X2) = \c2{xl + 2x1) + / Jc32:i(6x^ - xl) + ^C^{xl + 2xlf + V^. (2) 

The variables xi and X2 are connected to the deformation j3 and asymmetry 
7 by the standard relations [9] 

xi = /3cos7, a;2 = -^sin7. (3) 
v2 

As discussed in [8], the presence of bound states in atomic nuclei leads to 
a value of C4 > 0, whereas for C3 a positive value corresponds to a prolate 
shape, a negative value to an oblate shape. Similarly C2 may also be either 
positive or negative. 

The shape of the nuclear potential V{xi,X2) is a function of C2 and 
X = 0^/(0204). For C2 > 0, and < x < 56/9 the nucleus is spherical; for 
56/9 < X < 7 the nucleus is spherical in the ground state (g.s. spherical) and 
deformed in the excited states (e.s. deformed); for x > 7 it is g.s. deformed 
and e.s. spherical; for C2 < it is g.s. deformed and 7-unstable in the 
excited states. 

The chaotic behaviour of this model has been studied analytically by 
the authors of [9] using the criterion of negative curvature of the potential 

4 



energy. However, this criterion has been demonstrated to be inaccurate in 
[10]. The authors of [9] have also obtained the energy levels by quantizing the 
Hamiltonian (1) for g.s. spherical nuclei. For g.s. deformed nuclei, like rare- 
earth, however, it is not simple to obtain the correct energy levels because 
the potential energy is an asymmetric triple well. 

Table 1 shows the parameters of the nuclear potential (2), for the rare- 
earth ^^^Gd and ^^^Er according to [8]. The potential energy for these nuclei 
is plotted in Figure 1. 

3. The GFND method 

To obtain the energy levels, we use the Green function numerical diag- 
onahzation method (GFND) [11-14], whose starting point is the integral 
Schrodinger equation 

J J K{xi,X2, yi, y2] £)ipn{yi, y2)dyidy2 = e~^^"ijjn{xi, X2). (4) 

K{xi,X2,yi,y2',£) is the Euclidean short time propagator and ipn{xi,X2) are 
the Hamiltonian eigenf unctions. This equation can be approximated using a 
numerical integration rule which gives 

N N 
h=l k=l 

il^ij = ^Pn{xi,i, X2,j), Kfj^^ = K{xi^i, X2,j,xi^h, X2,k\ s), W2,i are the weights 
associated with the integration rule (a good choice is the trapezoidal one). 
The intervals of integration in Eq. (4) go from —00 to +00, but in Eq. 
(5) we take finite matrices and consequently finite intervals Li and L2. This 
corresponds to confining the system in a "box" of sides Li and L2. However, if 
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the intervals are large enough, this gives negligible corrections to the energies 
and to the wave functions of the bound states. Therefore, by diagonalizing 
the matrix Wi^hW2,k ^fjhk^ ^® obtain the energies and the wave functions 
directly. 

The Euclidean short time propagator can be written in the form (we use 
units B = h = 1): 

1 f 1 

K{xi,X2,yi,y2;£) = —exp l-—[{xi - yif + 2{x2 - y2f] - fixi,X2,yi,y2 

(6) 

where the first term in the exponential corresponds to the kinetic part of 
the Hamiltonian and the function f{xi, X2,yi,y2', s) is the potential term. Its 
explicit expression depends on the prescription chosen; for example the last 
point rule f{xi, X2, yi, y2', £) = eV{xi,X2) gives a propagator which is correct 
up to 0{e) only, while the symmetric expression - [V{xi,X2) + V{yi,y2)] is 
correct up to 0(£^). Moreover, a systematic expansion of the short time 
propagator in £, Ai = xi — y\ and A2 = ^2 — y2 is also possible [13]. This 
expansion gives rise to an asymptotic series, which, if used within the GFND 
method with non-singular potentials, allows one to obtain high numerical 
accuracy. There follows a brief sketch of the expansion method. 

Since the potential term of K{x\^ ^2, yi, ^2; is weighed by the Gaussian 
function: exp ~ + 2(^2 — y2)^]|, the difference Aj can be con- 

sidered of the order -y/i. As a consequence, we can expand in Ai and A2 
as well as in e. Therefore, following Ref. [13], let us write the short time 
propagator in the form 

1 1 

K{xi,x2,yi,y2;£) = ^^^Pi~^[(^i - 2/1)^ + 2(2:2 - 2/2)^] 



N 2{N-v) 2 2 

-EE E---E5;r''''(yi,y2)(x,, (7) 

v=G //=0 ii=l i^=l 

By expanding up to 0{e^) and requiring that X2,yi,y2', £) satisfies the 

Schrodinger equation, we obtain the coefficients g'^^l^'"''^'^ (see Appendix 1). 

The GFND gives accurate results and is more precise than techniques 
based on the direct diagonalization of the Hamiltonian, especially in the case 
of tunneling problems (for a deeper discussion of this point see [11-13]). 
Therefore, it is an appropriate method of studying the potential given in Eq. 
(2), where we have small tunneling effects between wells. 

4. Numerical results 

As is well known, the classical global instability of a system is appro- 
priately studied using Poincare sections. Regular regions are characterised 
by sets of invariant intersection points; chaotic regions by points which are 
distributed irregularly [15]. 

Prom the quantal point of view, the energy spectrum obtained using the 
GPND, has quasi-degenerate levels, due to the presence of three similar po- 
tential wells, two of which are identical. Therefore, it would be necessary 
to separate the energy levels into three classes, corresponding to the states 
localized in each well. Por our purpose, however, we need only take one of the 
two double degenerate levels, which are easily identified. Thus, we obtain 
approximately one hundred energy levels, which are enough for statistical 
significance. The energy spectrum thus obtained has been mapped into one 
with a quasi-uniform level density, by performing the unfolding procedure 
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described in detail in reference [16]. Then, the distribution P{s) of spacings 
between adjacent levels and the spectral rigidity A3(L) [17,18] has been calcu- 
lated. These spectral statistics are compared to Poisson statistics P{s) — 
and A^{L) — L/15 of integrable systems, and to Gaussian Orthogonal En- 
semble (GOE) statistics P{s) = fse"^"^/^ and A3(L) = n-^ln{L) - 0.0007 
of chaotic systems. 

Fig. 2 shows that in ^^°Gd, which has a high saddle energy (5.56 MeV), 
there is coexistence of regular and chaotic motion for energies below the 
saddle energy. In this region, the system has only a few energy levels and so 
we cannot perform a good statistical study of quantum levels. 

Fig. 3 and Fig. 4 show that in ^^°Gd above the saddle energy, at about 
8 MeV, there is chaos. For higher energies, at about 18 MeV, there is a 
quasi-regular behaviour. The non-universal behaviour of A^lL) for high 
values of L, not predicted by GOE, has been explained by Berry [19] using 
semiclassical quantization. 

Fig. 5 and Fig. 6 show that in ^^^Er, which has a low saddle energy 
(1.89 MeV), for energies above the saddle energy, at about 4 MeV, there is 
prevalently chaotic behaviour; for higher energies, at about 15 MeV, there is 
a predominance of regular classical trajectories. 

Conclusions 

In conclusion, the collective vibrational behaviour of the above nuclei 
displays a quasi-chaos quasi-order transition as a function of the energy. 
The mixed behaviour at very low energies cannot be easily shown by spectral 
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statistics because only a few energy levels are present. 

For the sake of completeness we observe that other phenomenological 
collective models have been recently proposed, see for example [20] and ref- 
erences therein. 
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Appendix 1 



The coefficients ' of Eq. (7) are given by (the yi and y2 dependence 
in g and V is understood) 



goi 



9l\ 



9T 



111213«4 

y4i 



Jll2«3M26 

i/51 



lll2«3l4»B«6 



5'02 
5'12 

9-n 



932 — 932 



932 — 932 



9ll' 



9l^ 



222 

932 



1 dV_ 

2 dxi^ 
1 



6 dxi^dxi^ 



24 dxj^-^dxj^Q^dxj^^ 



120 dxi^dxi^dxi^dxi 

1 a^i/ 



720 dxi^dxi^dxi^dxi^dxi^ 
1 (9^1/ 



5040 dxi-^dxi.^dxi.^dxi^dxi^dxiQ 

1. (V^(0,2) ^ ^(2,0)^) 

^ + ^(3,0)^ y{0,3) _^ y(2,l) 

^ \/(2,2) _^ y(4,0)^ y(l,3) _^ y(3,l) 
^ 1/(1,3) _^ y(3,l)^ y(0,4) _^ y(2,2) 

J_ /'T/(3,2) , t/(5,0)n 

360 ^ ^ 

^ ^y{2,3)^y(4,l)N 

360 ^ ' 

J_ (1/(1,4) ^ ^(3,2) ^ 

360 ^ ^ 
1 



1 

80 



360 
10 



(1/(0,5) _^ 1/(2,3)^ 



1111 ^ ^(04,2)^06,0)) 



2016 



2016 



„2211 _ ^2121 _ ^2112 

y42 — y42 — y42 

-g^-gr-gir - ^ (^^'"'^ + ^^'''^) 
g^ = g^'' = g^'' = g^' = ^ (v^''"' + v^'''^) 

gr - ^(^^°''^ + ^^'''^) 



5'03 = 



gls 



1_ (_0O,1)2^0O,4)/;lq 

24 ^ ' 

_ y(l,0)2 ^ y(2,2) /5 

1. (_v^{o,i)y(i,i) + y(i,4)/2o 
24 

+ 05,o)/2o + y(3.2)/io 

_ 01,0)^(2,0)) 
^ 1 (_0O,1)0O,2)^0O,5)/2O 

- i/(i.o)y(i.i)) 

^11 ^ _L (_01,l)2/9 _ 02,0)2/9 
_0O,l)02,l)/g 
_01,O)03,O)/8^02,4)/;^gg 
+ V(4,2)/84^y(6,0)/^g8) 
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+ y{i,5)/i68-y(i'i)F(2,o)/9 

-y(i'0)\/(2.i)/8 + y(3,3)/g4 

_p.(0,l)^{l,2)/g^^(5,l)/ig8) 

gll - ^(-^^^'^^V9 

_ 1/(0,2)2^9 _y (1,0) y(l,2)/g 
_y(0,l)y(0,3)/8^y(4,2)/;^g8 
+ 1/(2,4)/84^y(0,6)/;Lg8^ 

^04 = 1 (-1^(0-2)2/36 + y(6,o)/672 

_ V^(o,i)^(o,3)/;L2 + \/(°'6)/672 

_y(l,l)2/;^g_y(l,0)y(l,2)/;^2 

_y(2,o)2/3g_ ^(0,1) 1/(2,1) /;^2 
+ y(2,4)/224- v(1'0)1/(3,o)/;L2 

dx'^dy'^ 
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TABLE CAPTIONS 

Table 1: Numerical values of the parameters C2, C3, C4, Vq defining the 
nuclear potential V for the two rare-earth nuclei studied. 
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FIGURE CAPTIONS 



Figure 1: The potential energy of the rare-earth ^^^Gd (a) and ^^^Er (b). 

Figure 2: The Poincare section of ^^^Gd for the energy E — 4: MeV. 

Figure 3: Poincare sections for ^^°Gd with E — S MeV (below) and 
E^18 MeV (above). 

Figure 4: Spectral statistics for ^^°Gd with 6 < E < 11 MeV (below) and 
15 < E < 20 MeV (above); the solid lines are the Poisson statistics and the 
dashed hues are the GOE statistics. 

Figure 5: Poincare sections for ^^^Er with E — 4: MeV (below) and E — 15 
MeV (above). 

Figure 6: Spectral statistics for ^^^Er with 2 < E < 6 MeV (below) and 
13 < < 16 MeV (above); the solid lines are the Poisson statistics and the 
dashed hues are the GOE statistics. 
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Nucleus 


C2 (MeV) 


C3 (MeV) 


C4 (MeV) 


Vq (MeV) 




-100.64 


37.57 


668.01 


5.56 


166^^ 


-73.77 


64.96 


1138.52 


1.89 



Table 1 
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